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ABSTRACT 


This paper presents the modeling and simulation of wind energy Conversion 
System using the Permanent Magnet Synchronous Generator (PMSG). 
The objectives are: to extract the maximum power of the wind speed by 
controlling the electromagnetic torque of the PMSG, to maintain constant the 
DC-link voltage despite the wind speed variations and to attain the unity 


power factor. In order to ensure a regulation with high performance and a 
good robustness against the internal and the external disturbances, a new 
control strategy called the Active Disturbance Rejection Control (ADRC) is 
used. Therefore, the Analysis and simulation of the ADRC and PI controllers 
are developed with MATLAB/Simulink software. The performance of these 
controllers is compared in term of references tracking, robustness 
and grid faults. 
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1. INTRODUCTION 

In the recent years, wind energy became a famous clean energy source. The wind energy source 
systems are quickly growing thanks to the innovations and the researches in wind energy technology for 
electrical power production. Originally, a wind power source was used in standalone applications, however, 
wind systems connected directly to the grid [1]. Three types of generators are mainly usable in the wind 
energy conversion system: the squirrel cage induction generator (SCIG), doubly fed induction generator 
(DFIG) and Permanent Magnet Synchronous Generator (PMSG).For this work, the PMSG is used since it has 
small size, provides a self excitation and a high power density. In addition, the generator has the capacity to 
control reactive power, working also without the gearbox due to the higher number of poles for low speed 
operation [2, 3]. A configuration of a wind power system based on the PMSG is illustrated in Figure 1. 
The synchronous generator is connected to the grid through a back-to-back converter composed of the 
generator side converter and the grid-side converter [4]. 

The last years, many researches have appeared documenting the control of nonlinear systems related 
to wind energy. In the literature, proportional integral (PI) controllers based on field oriented control (FOC) 
is made up, this conventional technique presents some shortcomings, such as a low accuracy especially with 
the disturbance and nonlinearities according to dynamic characteristics of systems [5]. In order to overcome 
the disadvantages of the conventional PI controller which suffers from many limitations, a new control 
strategy called Linear Active Disruption Rejection Control (ADRC) founded by Han in 1995 is adopted to 
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remove disturbances affecting the facility [6,7]. This paper introduced an advanced method as the linear 
active disturbance rejection control (ADRC) based on vector control, in order to achieve a maximum power 
extraction, to maintain the DC bus at desired voltage value through the DC-link circuit, and to meet power 
stable of the wind energy conversion system [8]. 

In detail, the suggested design techniques propose a comparison between PI and ADRC technique 
which have been explored and used recently as a new technology for estimating and compensating 
uncertainties and disturbances, is based on Extended state Observer (ESO) which is the main part of 
the command. An evaluation is ensured to control the grid and generator side converter in order to regulate 
the electromagnetic torque and the stator current of the PMSG. Simulation is made under 
MATLAB/Simulink environment with a comparison between both commands. Results provided show clearly 
the feasibility and the effectiveness of the proposed two techniques. They present also a good characteristic, 
excellent output from the ADRC control comparing with PI command especially in terms of robustness and 
reference tracking. The sequential work flow of this paper is as follows: In section 2, complete architecture of 
the variable speed wind turbine conversion chain based on the PMSG has been described and the model has 
been determined. Section 3 covers both PI and ADRC design, followed by a control of grid side converter, 
bus voltage and generator side converter of the PMSG to maximize the production of the wind energy using 
respectively both controllers in section 4 and 5. Simulation results and discussion are made in section 6. 
Lastly, in section 7, a conclusion has been added to finalize the work. 
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Figure 1. Architecture of the variable speed wind turbine conversion chain based on the PMSG 


2. MODELING OF THE WIND ENERGY CONVERSION SYSTEM BASED ON THE PMSG 
2.1. Wind turbine 

The turbine transforms the wind energy into a mechanical energy available on a wind turbine rotor 
[9,10]. This power is expressed by: 


P, = =Cp(, B) pm R2V3 (1) 


Where p is the air density, R is the turbine radios and V is the wind velocity. The power coefficient 
Cp is given by the following expression [11, 12]: 


a5 
Cp (A, B) = Cy (52 — Cab — Ca) CF) 4 ca (2) 
where: 
A= 
1 1 0.035 (3) 


Ai A+0.08B 1483 


C\=0.5176, C2=1 16, C3=0.4, Ca=5, Cs=21, Co=0.0068 
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The maximum value of power coefficient Cp is Cpmax = 0.48, is achieved for pitch angle of rotor 
blades B is set to zero and Ag, =8.1, this point corresponds at the maximum power point tracking (MPPT). 
This model is translated by the following equation [21]: 


do. 
T, = Ty + fyO, +J- (4) 


Where Tg and Tem are respectively the mechanical and the electromagnetic torque of PMSG, J is 
the total moment of inertia and fv is the friction coefficient of the turbine. The electromagnetic torque 
reference is determined by the following equation: 


1 = 
Tem-ref = Fag, CPmaxP TR?! (5) 


The Figure 2 represents the block diagram of the turbine model with MPPT strategy: 





Figure 2. Model of the turbine with MPPT strategy 


2.2. Model of converters 
The generator side converter and grid side converter are controlled by the PWM control. These 
converters are represented by simple line voltages and the control signals Si (6) [8]. 


— 2S1—-(S2+tS3) 
Va = 3 Vac 


Í 
252-(S1+S 
t = ee Uae (6) 


| 
— 2S3-(S2+S1) 
(ve = SSD, 


2.3. Modeling of the PMSG 
The electrical voltages of the PMSG are given by the following equations: 


Vds = Rsias + La dt + €ds 
di (7) 
= Rigs + La -E + €e 
Vgs slags T Hq Gy qs 


Where the direct and the quadrature e.m.f components are expressed as follow: 
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| Cqs = —WsLgigs (8) 
eqs = WsLailas + Wye 


Where: 

La ,Lgare d and q axis inductances (H); 

Ris the stator resistance (Q); 

igssigs are the d and q axis machine current (A); 

dr is the flux linkage established by the permanent magnets in the stator windings (Wb); 

Ww, is the angular frequency of the stator voltage (rad/s); 

The expression of the electromagnetic torque is also expressed as a function of currents “(9)”: 


3 i 3 ioa 
Tem = 5 PP rigs + = P(La—La iasigs (9) 


p: Pairs of poles 
These equations have been grouped to form a simplified model of PMSG as depicted in Figure 3. 






TR: + SLa gs -3 / 2.p. By Ten 


-Vei : | 1/R; + SLq las” 
Figure 3. Simplified Theoretical Model of PMSG 


2.4. Modeling of RL grid filter and DC link 
The electrical voltages of the grid-side are given in the (d, q) reference frame by 
the following equations: 


digg 
& dt 
& dt 


Vind — R lq +L 
m gdg (10) 
Ving = Reglgg + L 

Where the direct and the quadrature e.m.f components are expressed as follow: 


(11) 


The active and reactive power supplied to the network can be expressed as : 


n = Vag ldg + Vag lag (12) 
Qg = Vag: lag — Vag: lgg 


The grid voltagevector is oriented on d-axis, then: 
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Vac = V, 

f dg E (13) 
Vag = 

Thus: 


Po = Vas.) 
l g dg: ldg (14) 


Qg = —Vag lag 


The dc-link capacitor is interface between generator side converter and grid side converter; By 
neglecting the converters losses, the state equation of the dc-link voltage are expressed as follow: 


Vac =$ J ic dt = =f (i, — i2)dt (15) 


where 
C: the DC link Capacitor. 


3. SYNTHESIS OF PI AND ADRC CONTROLLERS 
3.1. PI controller 

The proportional-Integral controller is a combination of both actions proportional and integral in 
order to cancel the static error. The structure of a parallel PI controller system is represented in Figure 4 
[10, 11, 16]. The derivative action of PID controller is excluded because it is characterized by amplifying the 
effect of system noise. 


ae 1/(Rs+s.L,) 
IqSret 





Figure 4. Structure of PI controller system for the PMSG 


3.2. Linear ADRC controller 
The mathematic expression of the first order for a time-varying dynamic system with single input u 
and single-output y is given as follow [13-16]. 


© = f(y,d,t) + bou (14) 
where: 
f(y, d, t) : the effect of internal and external disturbances, 
b0 : Represents parameter gain to estimate . 
d : The external disturbance. 
The mathematical system as follow describes the process [19-21]: 


X4 — X2 + bou 
| X2 — f (15) 
y = Xi 
f: = kX + Bi(y — £1) + bou (16) 
R2 = Boy — 1) 
Where the vector of the observer gain is expressed by: 
p= [sis 
= 17 
Bal — | S250 
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The pole of the observer SESO is determined by the technique of placement of the poles. Therefore, 
the control plant can be ensured by a simple proportional Kp (18). The input signal reference denoted r 


[20, 21]. 
uo = Kp- ĵ) = Kp (r — &) (18) 


where: 
Kp = SCL. 
SCL is the pole of desired closed loop. 
SESO = 3~7 SCL 
The only setting parameter is SCL. 


Figure 5 represent the Structure of ADRC controller and the ESO 





Figure 5. Structure of ADRC controller and the ESO 


4. CONTROL OF THE WIND SYSTEM BY PI CONTROLLER 
4.1. Control of the generator side converter 

The stator current is aligned on the q-axis (isd=0). Consequently, the electromagnetic torque is 
controlled by the q-axis current isq: The expression of the electromagnetic torque (9) become: 


3 : 
Tem = 5 PP figs (19) 


Thus: 
S (20) 


loa = —— hm 
qs-ref 2p, °™ ref 


The open loop transfer function is: 











s+ a 
O (kp +) gay) a (21) 


We use the method of poles compensation in order to eliminate the zero present on the transfer 
function [11, 12]. We pose: 


kj _ Rs 
nat (22) 
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Thus: 


kyo 
G(s) = ES e (23) 


Trs 


Where tr is response time: 


L 
w= (24) 


In closed loop we will have: 


1 
1+T;s 





F(s) = (25) 


Therefore, the parameters kp and ki of PI controller are given by: 


1 (26) 
k; = T Rs 


We chose: tr = 10 ms 


4.2. Control of the grid side converter and DC link voltage 
The open loop transfer function for the grid filter controlled by PI is: 





G(s) = (kp au )( : (27) 


Rg+tsLg 


Therefore, the parameters of PI controller are: 
(28) 


The deference between the active power generated by the wind generator Pm and the delivered 
power to grid P, is stored in the DC-link capacitor, the structure of DC link voltage control is shown in 
Figure 6. The converters losses are neglected, therefore the equation of the dc-link voltage is expressed by: 


Poc = Pm — Fg (29) 





Vd 
_ | 


Figure 6. Structure of DC link voltage control 


Figure 7 shows the PI control structure applied to the generator and grid side converter. 
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Figure 7. PI Control structure of the wind system based on PMSG 


5. CONTROL OF THE WIND SYSTEM BY ADRC CONTROLLER 
5.1. Control of the generator side converter 
The currents of the stator are rearranged as follow: 





Osa = = 1; + wo, 11, + SVa 
dt La La q La 
d e PS] _ La | a4 bay (30) 
dt lig sq Ws A sd La Pf ig sq 
These expressions can be written in the following form: 
Tsd = f(I,q,d,t) + bou(t) 
(31) 
Where: 
Rs Lq 
fisa d, t) = — L, sd + Ws La sq 
d d A (32) 
u(t)=Vsg and bọ=— 
La 
Where: 
Tsa = £(I,,,d,t) + boult) (33) 
dt ar y 
In the same way: 
Rs L 1 
F(Isq) d,t) = — L sq 7 Os E lsa ~ 7, OF 
q q ie (34) 
u(t) = Vsq and b y= T 
5.2. Control of the grid side converter and DC link voltage 
5.2.1. The currents of the grid filter 
The currents of the grid filter are rearranged as follow: 
dIag = 1 Rg 1 
w O dg Lg de e lag + [ma 
digg 1 Rey I 1 eo 
at -gp ag lag — 9s T lag + Ving 
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These expressions can be written in the following form: 


Zas = flag d, t) + boult) (36) 


dt 
Where: 


= Rg Lq 
fag. d, t) =m L, dg = Lg lag + Ws L4 lag 
i (37) 
u(t) = Ving and by =— 


E7 
Where: 

dI 

= = f(Iqg d, t) + bout) (38) 
In the same way: 


Za Rg La 
flag: d, t) ag Lp (98 = Le a8 Ts La 4g 
u(t) = Ving and bọ = es 


Lg 


(39) 
5.2.2. DC bus voltage control 
The power on the DC link Capacitor can be given by: 
dVac SEEE 
Pac = CVac a z Vac (i _ i2) (40) 


When the losses in the RL filter, the capacitor and in the two power converters are neglected, the 
powers across the DC bus is the difference between the generator power Pm and grid power Pg: 


Pac = Pm — Ps (41) 
Or: 

dVac 3 
CVac a = Vacli 7 5 Vaglag (42) 
Avago. Wes 3Vag . 
err rca (43) 


Taking A=Vdc?: 


dA 2vA. 3Vdg . 
eo ee (44) 
Thus: 
f(A, d, t) = yA l 
‘ 3Vdg (45) 
u(t)=igg and bo= ar 


Figure 8 shows the ADRC control structure applied to the generator and grid side converter. 
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Wind speed 


Generator Side Converter Grid Side Converter 








Figure 8. ADRC Control structure of the wind system based on PMSG 


6. SIMULATION RESULTS AND DISCUSSION 

As shown in Figure 9, the simulation was carried out with MATLAB/SIMULINK in order to 
validate the control strategy studied in this work. The system consists of several blocks, for example the 
turbine model with its wind profil, the model inverters with its parameters, and controllers based on the 
ADRC in Figure 10 and PI approach. All parameters are given in Appendix. A random wind profile is 
shown in Figure 11. To extract the maximum wind energy of the wind speed, an MPPT strategy is adopted. 
In order to maintain the DC-link voltage constant, a control of the grid side converter is necessary. This 
technique is based on the control of direct and quadratic grid currents, and ensures also an exchange of active 
and reactive powers between the stator of the PMSG and the grid. To attain the unity power factor, it’s 
necessary to regulate the grid reactive power to its desired value. In Figure 12, the rotor speed track the profil 
of the wind which allows the efficiency of the MPPT strategy. 








CLEESE E y 





Figure 9. The whole system simulation 
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Figure 10. ADRC block conception Figure 11. Profile of Wind speed 
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Figure 12. Rotor Speed of PMSG 


6.1. Test of reference tracking 

In Figure 15, the reference direct current Isa-ref is fixed at zero in order to get the unity power factor 
for the generator side converter (GSC). Figures 13 to 17 illustrate respectively the electromagnetic torque, the 
stator current, reactive power of grid and DC link voltage variations obtained by using the PI control and the 
linear ADRC controller. These parameters converge and tracks perfectly as their references (Tem-ref, Isd-ref, 
Isq-ref, Qg-ref and Vdc-ref), but with a remarkable static error and an important response time for the PI 
controllers. We can also notice that the response of these parameters is slow and presents some overshoots 
shown clearly in zoom figures. The main importance factor to analyze the performance of the ADRC 
controllers, according to fast variation of the references is described in Table 1. It shows a very good 
performance for the tracking test with high accuracy against to wind speed variation, fast dynamic response, 
a good stability comparing to PI methods, simplicity of design and implementation. 






Tem ref 


“O 5 


Figure 13. Electromagnetic torque of PMSG 
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Figure 14. q-axis stator current 
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Figure 15. d-axis stator current 
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Figure 16. Reactive power of grid 
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Figure.17 Vdc link voltage 
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Table 1. The system responses Parameters 
Type ofController PI ADRC 


Rise Time (s) 0.1 0.06 
Settling time(s) 0.15 0.07 
Overshoot (%) 1 0 


Steady state error (%) 5 0.02 


6.2. Test of robustness 

In order to investigate the robustness of the proposed ADRC control algorithm, two cases were 
selected in which we have changed the internal parameter of the PMSG. In the first one, an increase of the 
30% of the stator resistance nominal value is applied. Figure 18 and 19 shows results of a comparison 
between ADRC and PI controllers. It can be noticed that the characteristics have been regulated to its 
reference value. ADRC controllers present excellent performances with good efficiency, smaller overshoot in 
the reactive power of grid and DC bus voltage, fast response and neglected oscillation in the regulation of the 
direct stator current to zero. 





g 
i 
4 
=. = 
ü i *h 
4 
A 
=. m ; 
i 5 Tire i cect 1 i 
Figure 18. Electromagnetic torque of PMSG Figure 19. d-axis stator current 


In the second case, an increase of the stator inductance nominal value with 20% is applied. Figures 
20 and 21 illustrate respectively, the electromagnetic torque, the direct stator current. It is noted that the 
whole characteristics is regulated to its set value, and this results obtained by ADRC controllers are clearly 
more efficient than the classical PI. 


Tem ref 
Tem ‘Pr 
Tem ‘LADRC’ 





0 5 Time (sec) 10 15 





Figure 20. Electromagnetic torque of PMSG Figure 21. d-axis stator current 


6.3. Operation with grid faults (grid voltage dip) 

A symmetrical grid voltage dip is applied to evaluate the performances of our wind energy system 
controlled by ADRC strategy. Figure 22 shows that the voltage sag occurs at t=5s. DC bus voltage is 
controlled by ADRC and PI algorithm. As shown in Figure 23, the amplitude voltage oscillations is very 
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remarkable and present more than 1% comparing with the rated voltage in ADRC methods. If the voltage sag 
was very dip, a PI controller presents more oscillations which converge to unstable system. 
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Figure.22 Grid voltage 





Figure. 23 Va link voltage 


7. CONCLUSION 


In this work, a wind turbine conversion system using a PMSG connected to the grid were studied 
under two cases: normal and grid faults (voltage dip).To evaluate the performances of our system, a 
conventional PI and the active disturbance rejection control (ADRC) were compared. Simulations under 
Matlab/Simulink environment shows that the proposed ADRC methods gives very satisfactory characteristics 
with good efficiency, fast tracking and robustness for internal parametric variations of the PMSG, and 
external disturbances compared to the conventional PI controller even with the voltage dips effect. 


APPENDIX 


Table 2. PMSG Parameters [8] Table 3. Turbine Parameters [8] 


Rated power 6 kW Density of Air 1.225 kg.m3 
Pole pairsnumber p 5 Total moment of inertia Jt 0.41 kgm2 0.02 kg.m2 
Stator resistance R, 0.425 Q Optimal tip speed ratio Aopt 8.1 
Inductance La 0.0084 H Maximal power coefficient CPmax 0.48 
Inductance L, 0.0084 H Turbine diameter D 2m 
Rotor Flux ¢ 0.433wb 
Moment of inertia J 0.41 kg.m2 0.01197 kg.m? 
DC link voltage Uac 400 V 
Table-4. Grid and DC link Parameters [8] Table 5. ADRC Controller Parameters 
dc-link voltage 400 V Parameter gain by =1/L, 119.0476 
Capacitor of the dc-link 0.41 kg.m2 10mF controller gain K,= -Sct 400 
Grid frequency 50 Hz pole of ESO Sgso= 3ScL -1200 
AET ETEA Bo O Extended state observer gains (ESO) ne : ee 


Table 6. PI Controller Parameters 
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Parameter gain bo =1/L, 119.0476 

controller gain K,= -Sct 400 

pole of ESO SEso= 3Scr -1200 
B1=2400 


Extended state observer gains (ESO) B.=1440000 
= 
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